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Abstract 

In this paper we present exactly solved extended quintessence models; 
furthermore, through a dynamical effective Q-cosmological "constant", 
we recover some of the A decaying cases found in the literature. Finally 
we introduce a sort of complementarity between the Q-dominated or A- 
dominated expansions of the Universe. 



e-mail addresses: 
dcritis@na.infn.it 
marino@na.astro.it 



2 



There are recent crucial informations coming from obeservations which seem to 
put new light on our nowadays knowledge of the Universe. We mention just two 
of them: first, the recent observations on Type la supernovae are strongly in 
favour of an accelerated Universe [EJ Q . As well, measurements of microwave 
background, mass power spectrum M and lensing statistics (for example, [5]) 
all suggest that a large amount of energy density in the Universe should have a 
negative pressure. That is, we have to include some missing part of the energy 
needed to reach the critical one. Second, as suggested by BOOMERanG ||, 
it seems to be very reasonable to assume that the Universe is spatially flat, 
that is f2fc = (in agreement with inflationary scenarios) as we will consider 
all along the paper. In connection with the first information we will follow the 
assumption present in the literature assuming that the missing part of the en- 
ergy density can be supposed to have the form of "quintessence" or the form 
of a cosmological constant Q || |J [fio| . Actually, to both these possibilities 
is connected a negative pressure, in order to explain the above quoted observed 
acceleration: then we have to require such a component to have a state equation 
p = wp, with —1 < w < (very recent considerations Q require — 1 < w < —0.6 
for quintessence, see also |0, Q], Q, for the a;- matter scenario); when 
w = — 1 we have the cosmological constant scenario Jll| . 

In this letter we clarify some questions connected with these different models 
improving the results already published In the context of non- minimally 

coupled quintessence theories (quintessence is, from this point of view, directly 
related with geometry (l6)), we give an exact treatment of the models pre- 
senting one of the most commonly used quintessence potential, i.e. the inverse 
power potential V = VoQ~ a , with a > (for example see JlTj Q). We, also, 
solve exactly the model connected, in the minimal coupling case, to another 
type of potential now present in the literature, i.e. V(Q) — Vbsinh(aQ), p0[ . 
In this scenario we shall discuss a definition (already given in a different con- 
text [^l]) of quintessence effective cosmological constant, and then, using it in 
the cases we are considering, we explain some of the most ad hoc A decaying 
behaviours considered in literature [p6| . We will spend, also, some words on 
the (A-)fine tuning jio) and the cosmic coincidence problems (fl m ~ ft\) (see 
for example Jl5|]), which are present in these scenarios. We believe that such 
an exact treatment of the models can put some further light on the two above 
quoted problems connected with these scenarios. Furthermore we want to men- 
tion, in these introductory remarks, that, concerning the use of A to explain 
the supernovae results, doubts have been presented: for example it has been 
shown that the inhomogeneity can be used to explain them Together 
with this, Riess et al. [^3| have opened some doubts on the interpretation of 
same data stressing the possibility of interpreting them entirely in the context 
of a standard open Friedmann-Robertson- Walker model joined with a reason- 
able astrophysical evolution model of the white dwarf supernovae progenitors. 
Finally we will exhibit a sort of complementarity relation showing the usual 
mutually exclusive presence, in these models, of the quintessence evolving dy- 
namics and of the given (effective quintessence) cosmological constant evolving 
dynamics. This result, together with finding a very well studied A-decaying 
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from our quintessence denned effective cosmological constant, we believe is the 
bridge which makes evident the difference in using a quintessence dominated 
Universe or a quintessence defined effective A dominated Universe. This means 
that the last one can be " deduced" from the former one and, in a dynamical way, 
it could be connected to an asymptotic really constant A at the end (this last 
aspect is interesting also because we think it can reintroduce a new attention in 
the cosmological NoHair theorem p7j ) . 

We start considering the general (field) action that describes the model we 
are going to use: 

A = J V=g(F(Q)R + \grQ,„,Q, v - V(Q) + C m )d 4 x (I ) 

being (F(Q),V(Q)) the (generic) functions describing respectively the coupling 
and the potential, R is the curvature scalar, —g^Q^Q.v is the kinetic energy of 

the quintessence field and C m describes the standard matter content. In units 
8irG = ft, = c = 1 we recover the standard gravity when F is equal to -5. In the 
flat Friedman-Robertson-Walker cosmologies, (I) gives rise to the "point- like" 
Lagrangian C (we use this expression intending that the field density Lagrangian 
connected to (I), because of the cosmological principle, can be considered as 
defined in the minisuperspace where the remnant two field variables (a, Q) have 
to be considered functions only of the cosmological time and then considered as 
describing a two degrees of freedom model) . In this sense C is defined on a two 
dimensional configuration space (actually we are treating the quintessence field 
which is not homogeneous, as a function of the time only) 

C = 6Faa 2 + QF'Qa 2 a + a 3 p Q - Da~ 3h ~ 1) (2) 

where a is the scale factor, pq = —Q 2 — V(Q) and 7 is given by using the 

standard matter p m = (7 — l)p m state equation (completely independent of 
the second order differential Eqs. connected to (2)). Prime denotes the (total) 
derivative with respect to Q, dot the same with respect to time. 
From Lagrangian (2) we get the same equations which come from the field Eqs. 
derived from (1) in the FRW metric. Variation of Q gives the Klein-Gordon 
equation, whereas the Bianchi identity gives rise to (standard) p m = Z?a -37 , 
being the constant D given by D — /9 TOo Oq 7 (the expressions (.)o denote the 
value of the quantity (.) now). We find from (2) 

H 2 + -H+^ + ^= 0, (3) 
F 6F 6F ' y ' 

a 9 F F 1 1 
2 a +H + F +2H F-2F P *-2F P - = °> & 

where pq = —Q 2 + V(Q). In this presentation the Klein-Gordon equation is 
derived from (3), (4); it is interesting that Eq.(3) is precisely Ec = 0, being Ec 
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the (constant) energy associated with C. We will consider here only the dust 
case, i.e. p m = 0. According to Noether theorem (the Noether symmetry is 
studied in the quintessence minisuperspace) , we could get a further information 
(the possible existence of that symmetry) very useful to exactly integrate the 
system (3), (4) as well as to find a form for the two unknowns (F(Q), V(Q)). 
The study of the existence of this symmetry actually leads to an infinite set of 
Lagrangians (2) if the following relation between the two functions F(Q), V(Q) 
is satisfied: 

V = V (F(Q)) 2 ^ S \ (5) 

where p = — — ^ , and if F(Q) is of the form 
2s + 3 ' 

F = F (s)Q 2 , (6) 

with F (s) = 4S ( s +^f s+2 ) ■ The parameter s labels each Lagrangian belonging 
to the class of infinite Lagrangians (of type (2)) admitting a Noether symmetry. 
The value s — is a permitted value and gives rise to a model admitting a 
Noether symmetry, but has to be treated in a different way (for a complete 
treatment of this approach see Q). Requiring F{Q) < (attractive gravity) 
and V(Q) of inverse power-law type, we get s e (—3/2, —1). In Fig. (1) we 
plot the two functions F (s), p(s): (pratically all the interesting expressions of 




Figure 1: In (a) we plot the coefficient of the function F(Q). In (b) we show all the 
possible exponents for the inverse power-law potential: pratically, all exponent 
values are avalaible. 

the inverse power law potential are available: for example, V = is relative 
to the model we pick up fixing the value s = —1.257571). The existence of 
the Noether symmetry gives a further first integral of the second order Euler- 
Lagrangian Eqs. related to (2). Using it we can exactly solve those equations; 
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their solutions a(t), Q{t) are 

s + l 

a(t) = 6 2 (s)[kit + k 2 }£* {ft 3 [fcit + *2]^ + M + &o} 2 (7) 
Q(t) = Siis) j[M + fc 2 ]^ {fc 3 [fcif + fc 2 ]^ +6^ + 60! 3 1 (8) 



being 

5x( S ) 



X( S ) X ' 



H^7 



, -«D . s + 3 S 42 Vb(s + 3) 



7(s) 



3E ' a 7(«)' 3fcfr(s)(« + 6) 

2s + 3 
12(« + !)(« + 2)' 



where So, wo, &o and 61 are the four initial data (the constant So conies from the 
existence of the Noether symmetry); considering the condition Ec = 0, we get 
3So&i + sD = 0, representing the only constraint on the four initial data for the 
system of the two second order differential Eqs. (7), (8) which, then, become 
three initial data as usual. From this constraint we see that neither S nor 
61 can be zero because we are studying models with nonzero standard matter. 
The constant D comes from the Bianchi identity for the standard matter as 
we have already stressed, and it cannot be considered like an initial datum for 
the system (7) and (8) because the state equation is used. These informations, 
together with Vq ^ 0, which is quite obvious, tell us that k\ and £3 have to be 
different from zero. It is important to stress that > 1, for s G (—3/2 , — 1); 
then, for large t the two functions (7), (8) become a(t) = A a t r{ - S \ Q = Q a t N(s \ 
where 

fi -I- Qs -I- 9 s 2 

rW = T(3 + s) (> °' f ° r S G M/2 ' ' (9) 

N(s) = ~( 3 + 2s ) ( >0 , for se (-3/2,-1)), (10) 
s 

which we plot in Fig. 2. We have indicated with (A ,Q ) the two coefficients 
which come from (7), (8) for f > 0: these two constants are parametrized by 
s, that is they are dependent on the model; they also depend on S and k\. It 
is important to stress here that this does not give rise to any real limitation. 
Actually, we are just facing a typical situation found in standard cosmology. 
If we consider, for example, the de Sitter solution which is given by a(t) = 

aoe"^ 2A / 3t , of course we have to impose ao > from the very beginning, and 
this does not imply that such a behaviour depends on initial data. Actually this 
is required by considering an expanding Universe; then the initial data have to 
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Figure 2: The plots of the two exponents we find in the large t behaviour of the two 
functions Q and a respectively. It is interesting to see the these two functions give 
rise to a monotone t-dependence for both the functions we are considering: then the 
inverse functions (say t = t(a)) are always well defined. 



belong to the semiplane ao > 0. The same concerns (Ao,Qo), that is for any 
of those data the asymptotic behaviour of a(t) and Q(t) is the same, they both 
cannot be zero for all the physical initial data one can choose. We believe that, 
in this way, we do not have any (initial data) fine tuning in our models, this 
being strictly connected to the control we have on those data; that's why is 
so important to have the exact solutions of the evolution of the cosmological 
variables. It is also relevant that we find most of the (large t) solutions spread 
out in the literature. We see that Q(t) diverges for any s S (—3/2, —1), the 
same is true for a(t): in the same s-interval Q goes to zero for large t. 
Of course, we can fix a special value of s (which identifies the model) in order 
to get special time dependences for the scale factor: for example, we can have 

a ~ i 2 / 3 (11) 

% 

for s = —1.0788. For this value of s, the potential V(Q) is given by 56 , 
(Vo = F (s) 2 p^Vo). 

Q(t) diverges like i - 790 (for s = —1.0788). In Tablel are given particular values 
of the s-functions which play a relevant role (actually, we recover almost all of 

the studied beaviours for a(t)). Furthermore we see that, being fin — —H + ^Q, 

1 

and substituting the solutions we found, we have on ~ -rr = „ , , in the case 

s = —1.0778- model, we get r = 3/2, and on ~ — j. Then, among the infinite 
models under considerations, there is an exactly integrated one which shows a 
scaling-type behaviour j24[], fin ~ — scales like p rn (we also see that there are 



5 



s 


r(s) 


l/r(s) 


N(s) 


1/N(s) 


2p{s) 


-1.0788 


0.666 


1.5001 


-0.780 


1.280 


-0.280 


-1.1 


0.708 


1.412 


-0.727 


-1.375 


-0.75 


-1.2 


0.888 


1.125 


-0.5 


2 


-2. 


-1.3 


1.0497 


0.952 


-0.307 


3.25 


-4.5 


-1.4 


1.196 


0.8358 


-0.142 


7 


-0.280 


-1.499 


1.332 


0.750 


0.001 


749 


-4.5 



Table 1: For different choices of the parameter s we give the values of the functions 
describing the most important behaviours appearing in each models we consider. It is 
interesting to stress that for s = —1.499 we get the radiation behaviour for the scale 
factor. 

a large number of models exhibiting 2/r < 3 (tracker behaviour, see ]Tot). 
In Tablel we find also the relative values for the power in the expression of the 
potential. 

The treatment reported before requires s ^ |l]|] ; as we said the case connected 
to the value s = has to be discussed separately. It can be shown that two 
subcases are found: there is a Noether symmetry for the Lagrangian (2) if: 

(i)F = K Q 2 , (K < 0) and V = V Q Q 2 (V > 0); (ii) F = -- (minimal 

coupling), and V = V {Ae» Q - Se^) 2 , with fx = ^/3/2 |(J. 

It is easy to see that in this last case the starting Lagrangian we have to con- 
sider is C — 3aa 2 — a 3 (^Q 2 — V{Q)) — D for the standar dust case. Everything 
goes like the s ^ 0-cases and we will not discuss here the problems related to 
the initial data as well as the constraint ( Einstein (0, 0) Eq.) imposes on the 
four initial data. We will discuss only the minimal coupling case which splits in 
two subcases (for the complete discussion of both cases, see Jlj]]). Depending 
on the relative signs of the two free parameters (A, B), the general solutions of 
FRW Eqs. are: 
AB > 



o(i) = 

Q(t) 



w (\/ABt + z Q ) 2 - u 2 sm 2 (^/ABt + wi) 



3 V A 



4AB 
{uoVABt- 



1/3 



zo) 



■ u)q sin(V ABt + uji ) 



(woV ABt + z Q ) — ujq sin(V 'ABt + uji) 



(12) 
(13) 
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AB < 



Q(t) 




(w y/-ABt + z ) 2 + ujI sin h 2 (^-ABt + u>i) 
i\AB\ 



1/3 



(uj y/-ABt + z ) + ujQ sin h(y/-ABt + uj 1 ) 



-(uoV-ABt + z ) + ujq sin h(y/-ABt + oji) 



(14) 
(15) 



where u>q, lu\, zq are the three integration constants (of course A, B and D are 
different from zero). In both cases, Q — > const, for large t; in the first case 
a(t) behaves as t 4 / 3 for small t and as i 2 / 3 for large t (self tuning solution sec 
[p4|), in the second a(t) has a de Sitter asymptotic behaviour. Anyway, the 
asymptotic behaviours are independent of initial data in the same sense that we 
have clarified above. 

After having exhibited exact solutions for inverse power law potential, and for 
some kind of exponential potential, let us go now to discuss a way of introducing 
an effective, time dependent, cosmological "constant" [^lj (we hope in this way 
to solve the connected fine tuning problem j30|]). Before presenting our defini- 
tion of an effective, cosmological "constant" it is notewhorty to recall that: i) 
standard A is introduced by hands; ii) it determines the (in general asymptotic) 
time behaviour of a(t) through the (0, 0)-FRW Eqs., which can be rewritten as: 

(h — (^H + = g m . An expanding universe requires, asymptoti- 

cally, H = \ — . If we look at our (3) we see that it is still possible (in the dust 



case) to put it in the similar form (H — Ag^ ) {H 
being 



l^effj)--^ 



Peff' 



A, 



effi 



F 
2F 



\ 




PQ_ 
6F 



(> 0, because F < 0) . 



(16) 



The second root, i.e. A e g- 2 , is less than zero and does not affect the asymptotic 
time behaviour of the cosmological quantities. Two comments on definition 
(17) are in order: first, it is completely defined on the quintessence side of the 
quintessence-tensor theories we are considering, i.e. it is defined only using Q 
once we have F(Q) and V(Q); from this point of view, A e g- is not introduced 
by hands but using the same procedure we have mentioned above, that is using 
the roots of the (0, 0) Einstein Eqs. We, also, want to stress that it could be 
of some interest to study, in the effective cosmological constant scenario, under 
which conditions (17) becomes, for t>0, a constant giving back the condition 
for having a de Sitter asymptotic behaviour (the complete discussion of how the 
asymptotic cosmological NoHair Theorem can be generalized to this case is in 
pl|). From definition (17) we have that, using the solutions of the system (3) 
and (4), Ag^ is a function of (Q(t), Q(t)), and of the parameters connected 
with F(Q) and V(Q). It is noteworthy that for s 6 (— 3/2,,— 1) it is found 
that Agg^ decays in a way very well studied in literature, even if in all these 
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discussions the decaying of A e g- is given ad hoc (see the exhaustive paper by 

Overduin and Cooperstock pq|). 

Using the solutions (7), (8) for t » 0, we find: 

B 1 

iYpff, = — m = m(s) = — — . (17) 
eiii am w r ^ v > 

Even if the constant B is depending on s and on the initial data of the problem 
as well as on the parameter Vq, as we have stressed till now, the way 
behaves in time does not depend on the initial data. Then Eq. (18) shows that 
the form of the decaying is general, that is A e g- i goes to zero independently of 
any initial data. In |p6| the form (18) is given a-priori and then an explanation 
of the constant B is not found (in terms of H we have also that A e g ~ H). 

When a(t) ~ t 2 l z we have m — 1.50016, that is Agg^ ~ compatible with 

lens statistics, power spectrum of matter density perturbations f26| . Also the 
value m = 2 is permitted for s = —1.0002. In the case s = 0, with AB > 
B 

we have that Agg^ ~ 3 ^ 2 , [B is the equivalent of B introduced above, in the 

case s = 0, Q goes to const.). For AB < we get the interesting result 
Agg^ —* A(A, B) — const., independently of any initial data and then we recover 
a de Sitter asymptotic behaviour for the scale factor (also in this case Q goes 
to const.). 

In summary, in the context of Noether symmetry approach to non-minimally 
and minimally coupled quintessence tensor theories of gravity we have obtained 
a class of models exhibiting two important types of potentials, respectively, 
inverse power law and exponential, and we have exactly solved it. If we look at 
the solutions, we can state a property which shows a sort of complementarity 
between quintessence and the dynamically defined A e g. Before writing down 
this property, in Table2 we give the behaviours we have found for (a(t), Q(t)), 
Q(t)) and for the effective A-term we have introduced. The relation we propose 
is the following: 

A effl /(Q))=cost., t»0, (18) 

where / depends not only on the Q-field but is determined also by s (for exam- 
i 

pie is given by Q N ( s ) in the cases s € (—3/2, —1), i.e. it is determined by the 
model we choose. Relation (19) is true, in general, if we assume the esistence 

of Noether symmetries and that |J = G e g-/G e g- ^5 (see Jl9| for a better 

understanding of the hypothesis on -p = G e ff/G e g 0). We see that (19) 
gives rise to the already mentioned complementarity between the A-term and 
the (divergent or convergent to a constant) quintessence content in the universe. 
From this point of view we do not have that dark energy, or quintessence, and 
cosmological constant are completely different forms of energy: first, because 
they are different forms of the same (quintessence) energy; in fact, we derive the 
A-term from the presence of quintessence. Furthermore, they are complemen- 
tary in the sense given by (19) (actually, the quintessence kinetic energy has no 
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s e (-3/2, - 1) 


s = U, An > U 


s = L), An < 


a(t) ~ 


a(t) ~ i 2 / a 


a(t) ~ exp[A effi i] 


Q(i) - oo 


- const.(^= 0) 


Q(t) ~ const.y^ 0) 


Q(i) - 


Q(t) - 


Q(t) - 


F(Q(t)) ~ oo 


F(Q(t)) = -1/2 


F(Q(t)) =-1/2 




F/F = 


F/F = 


V(Q(t))~0 


V(Q(t)) ~ 


- const.(^ 0) 


A effi ~ 1/a™ 


A cff 1 ^ A 


A e g ~ const. (7^ 0) 



Table 2: The t > behaviour of all the physical quantities introduced. It is notewor- 
thy that in all the cases the potential rolls down to its minimum. In the AB > case 
its minimum is zero, whereas in the AB > case the minimum is different from zero, 
and then we recover a real cosmological constant. We see also that the kinetical-Q 
energy has no asymptotical role 

asymptotic role in all the cases we have presented) , from which we deduce that 
the dominant ingredient is connected with the large t behaviour of the Q-field 
and with the form of the potential minimum. When this minimum is zero we 
have that the Q-field dominates, when this minimum is non zero we have the 
asymptotic dominance of the effective cosmological " constant" . In this paper 
we do not discuss the way A, as well the solutions we have found, depends on 
the potential parameter; to this purpose see our Jll|. From Table2 it is clear 
that the two possible evolutions of the scale factor are quite different, then the 
expansion history of the Universe is different: its acceleration and its age will 
strongly depend on the dominance of the Q-matter or of the effective A-Q. In 
all cases the potential rolls down to its minimum and then, as we have stressed, 
the value of this minimum plays a very important role: in the cases we have 
presented, when this minimunm is different from zero, we have an asymptotic 
cosmological constant (de Sitter behaviour for the scale factor) ; in the cases this 
minimum is zero, we do not recover any cosmological constant (power expansion 
of the scale factor). More precisely, we have that in the (Noether) nonminimal 
coupling case the minimum is always zero, whereas the coupling is divergent; 
then we get that G e fj/G e g- and that A decays. In the (Noether) mini- 
mal coupling case, when AB > 0, the potential minimum is zero, whereas in 
the AB < case the minimum is different from zero, and then we recover an 
asymptotical true cosmological constant. In fjq| we discussed in more details 
the possible values of the parameter w, finding a range of values, at least for the 
case s £ (—3/2, —1), coherent with the currently used range of values as given 
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in the introductory remarks. We have to say that, concerning w, compared to 
what is found in the literature, we have followed here a different approach: we 
have not introduced this new parameter but we have reconstructed it using the 
knowledege of the exact solutions: more precisely the state equation pq = pq{q) 
has been found knowing pq = (t(pq)), that is using a(t) and Q(t) for t ^> 0. 
Concerning the large scale structure connected with this approach, see our pa- 
per pl[ . We believe that, through the knowledge of the exact solutions we have 
presented, it is possible to have a better control of the role of initial data; ac- 
tually, all along the paper we have shown that using our approach we do not 
have that problem. Furthermore, using our approach, we hope to have clarified, 
among the reported other features, some aspects of quintessence models: more 
precisely we have given relation (5), we believe interesting, between the coupling 
and the potential as well as the complementarity relation (19) which holds for 
the types of couplings and potentials we have studied. 
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